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We report on a calculation of the light hadron spectrum and quark masses in three-flavor dynamical QCD 
using the non-perturbatively 0(a)-improved Wilson quark action and a renormalization-group improved gauge 
action. Simulations are carried out on a 16 3 x 32 lattice at /3 = 1.9, where a -1 ~ 2 GeV, with 6 ud quark masses 
corresponding to m„-/m p ~ 0.64-0.77 and 2 s quark masses close to the physical value. We observe that the 
inclusion of dynamical strange quark brings the lattice QCD meson spectrum to good agreement with experiment. 
Dynamical strange quarks also lead to a reduction of the uds quark masses by about 15%. 



1. Introduction 

As a joint collaboration of CP-PACS and 
JLQCD, we have initiated a project for the pre- 
cise measurement of the light hadron spectrum 
and the quark mass with three-flavor full QCD. 
The systematic deviation of the quenched lattice 
spectrum from experiment found in is signif- 
icantly reduced with the inclusion of dynamical 
up and down quarks |2I3| . We now extend the 
work to three flavors. 

We adopt the renormalization-group improved 
gauge action and the 0(a)-improved Wilson 
quark action with the non-perturbativcly im- 
proved clover coefficient csw obtained in jlj . Sim- 
ulations are carried out at /3 = 1.9, at the scale 
a" 1 = 2.05(5) GeV, on a 16 3 x 32 lattice us- 
ing a Polynomial Hybrid Monte Carlo algorithm 
developed in 0- Six values of the ud quark 
mass, corresponding to tops,ll/"1v,ll — 0.64- 
0.77 1 , are used to extrapolate to the physical 



value m^/irip ~ 0.18. Two values are taken for the 
s quark mass, corresponding to mps,ss/ m v.ss — 
0.72 and 0.77. They are close to the semi- 
empirical value m 7?s /m0 = 0.68 from chiral per- 
turbation theory, so that a short extrapolation is 
enough to get to the physical value. 3000 HMC 
trajectories are accumulated at each simulation 
parameter. 

2. meson masses 

We focus on the meson spectrum since our lat- 
tice size of La ~ 1.6 fm is too small for baryons 



Fig. ^ shows the J parameter |B] defined by 
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x In this report, subscripts LL, LS and SS represent the 
light-light, light-strange and strange-strange mesons, re- 



This is compared with the previous results in two- 
flavor and quenched QCD [2], which are evaluated 
for the physical kaons mpg/my = wikIvtik* a s 
J = mv(dwv/dmp g ) after taking the chiral limit 
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Figure 1. The J parameter in two- and three- 
flavor QCD (left panel) and in quenched QCD 
(right panel). 

for sea quarks 2 . We see that J in three-flavor 
QCD stays at values systematically higher than in 
quenched and two-flavor QCD, and is consistent 
with experiment. 

Our lattice size La ~ 1.6 fm is significantly 
smaller than that used in quenched or two-flavor 
QCD. Because finite-size effects generally lower 
the J parameter we conclude that the consis- 
tency of J between three-flavor QCD and exper- 
iment is not an artifact of finite-size effects. 

Chiral extrapolations are carried out using 
polynomial forms. We fix the lattice spacing and 
the ud quark mass with m w and m p . The strange 
quark mass is fixed with either ttik (X-input) or 
?7i0 (0-input). 

In Fig. |3 is plotted as a function of the 

lattice spacing. The three-flavor result is higher 
than that in two-flavor QCD at comparable lat- 
tice spacings. The consistency of the three-flavor 
value with experiment seen already at a finite 
lattice spacing lends us hope that scaling viola- 
tion with the use of the non-perturbatively 0{a)- 
improved action is indeed small and the consis- 
tency with experiment holds in the continuum 
limit. This should be verified in a future study. 

Similar consistency with experiment can also 
be seen in other meson masses as shown in Fig. [3J 
where deviations from experiment are plotted. 
Improvement over the quenched results of 1 is 



2 We now prefer the definition of J by Eq. (T\ . as it is 
more straightforward, using the hadron masses measured 
in the simulation without resorting to any extrapolation 
procedures. 
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Figure 2. K* meson mass as a function of lattice 
spacing. 
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Figure 3. Deviation of meson spectrum from 
experiment in three-flavor (filled symbols) and 
quenched QCD (open symbols). 

evident, as anticipated from the agreement of the 
J with experiment. 

3. Quark masses 

The physical quark masses, m u d and m s , are 
calculated through the axial- vector Ward identity 
(AWI): m q = (A 4 AiP)/(2 (PP)). The matching 
to the MS scheme is carried out at /_i = a~ 1 using 
the one-loop Z factor with the tadpole im- 
provement. The renormalized mass is evolved to 
/i = 2 GeV using the 4- loop j3 function. 

The vector Ward identity (VWI) may offer an 
independent determination of m uc i and m s . The 
critical hopping parameter K c for ud quarks, de- 
fined where the light-light PS meson mass van- 
ishes, depends on m s due to the lack of chiral 
symmetry. The two natural choices, K c (m s = 
0), i.e., K c in degenerate three-flavor QCD, and 



K c (m s — mf ys ) with the physical m s , however, 
lead to m u d with the sign opposite to each other, 
indicating a large uncertainty from chiral symme- 
try breaking at a -1 ~2 GeV. The determination 
of the VWI m u d requires a careful study of scal- 
ing violation; we focus on the AWI quark mass in 
this report. 

Fig. 01 compares m U d an d m s in three-flavor 
QCD with our previous results in two-flavor and 
quenched QCD |1I2| . The discrepancy in m s be- 
tween K- and 0-inputs in quenched QCD disap- 
pears in two- and three-flavor QCD. Note that 
m u d and m s decrease with an increasing number 
of flavors of dynamical quarks. 

Taking the AWI mass with the if -input as the 
central value and including the deviation between 
the K- and 0-inputs into the error, we obtain 
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in three-flavor QCD. Dynamical strange quarks 
decrease both m u d and m s by about 15%. The 
ratio m s /m uc i is consistent with the one-loop es- 
timate of chiral perturbation theory 24.4(1.5) [HJ- 



Figure 4. Comparison of m u d (top figure) and 
m s (bottom figure) between full and quenched 
QCD. The left, middle and right panels in each 
figure show results in three- and two-flavor and 
quenched QCD, respectively. 



4. Conclusion 

Our simulation for three-flavor QCD shows 
that the inclusion of dynamical strange quarks 
brings the lattice QCD meson spectrum to good 
agreement with experiment. A non-negligible re- 
duction is also found for quark masses. These 
preliminary results are sufficiently encouraging 
to warrant further detailed and careful investiga- 
tions with time-consuming three-flavor full QCD. 
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